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The greater confinement of light that is possible in photonic crystal fibres leads to a greater 
effective nonlinearity, which promises to yield greater quantum squeezing than is possible in standard 
optical fibre. However, experimental work to date has not achieved improvements over standard 
fibre. We present a comprehensive numerical investigation of polarisation squeezing in photonic 
crystal fibre in a Sagnac configuration. By including loss, a non-instantaneous Raman response, 
excess phase-noise, second- and third-order dispersion and self-steepening, the simulations are able 
to identify the physical factors that limit current photonic crystal fibre squeezing experiments. 



Squeezed light is an important resource for the fields of 
quantum information [1, 2], where it is a route to quan- 
tum entanglement [3-5], and precision measurement, 
where for example it is being used to increase the pre- 
cision of gravity wave detectors [6]. It can be generated 
simply by an intense light pulse propagating in an optical 
fibre [7], via the optical Kerr effect [8] which is induced 
by the intrinsic third-order or nonlinearity. 

The advantages of fibre based squeezing include its 
flexibility of wavelength and power, avoidance of the 
phase-matching required in x^ materials, low intrinsic 
losses, relatively low cost, long form factors and ability 
to integrate with existing optical network infrastructure. 
It also provides a source of continuous variable entangle- 
ment [2], useful in unconditional quantum information 
experiments. The biggest disadvantage is two dominant 
noise sources that inhibit the degree of squeezing. Raman 
noise scales inversely with pulse width and guided acous- 
tic wave Brillouin scattering (GAWBS) scales with fibre 
length. The pursuit of better squeezing in optical fibre 
becomes an exercise in controlling these noise sources. 

Polarisation squeezing has to date achieved the best 
result in optical fibre, with 6.8 dB [9] below shot noise. 
One promising solution to curtail the noise sources is 
photonic crystal fibres (PCF). The higher confinement 
in PCF means that the same Kerr effect can be achieved 
with shorter interaction lengths or lower pulse intensity 
than standard optical fibre. A shorter fibre length means 
a shorter distance for classical GAWBS noise to accu- 
mulate. On the other hand, a less intense pulse, with 
a broader pulse shape, will not excite the Raman re- 
sponse to the same extent. Despite this potential, im- 
proved squeezing with PCF was been difficult to achieve 
experimentally [10], with 3.9 dB below shot noise the best 
achieved to date [11]. 

In this paper we present the results of comprehensive 
numerical simulations of the PCF squeezing experiments 
of Milanovic et al [11]. The physical factors that pre- 
vented the PCF experiments from achieving the same 
levels of squeezing as those achieved with standard fibre 



are identified to be an excess classical phase noise and 
an additional effective loss arising from reduced fringe 
visibility. 

The experiment employed a Sagnac configuration of 
counter-propagating pairs of pulses, which was found to 
be an improvement over the earlier single-pass polarisa- 
tion squeezing [10]. After a 50:50 splitting, ultrashort 
soliton pulses counter-propagate through a lm length 
PCF in the same polarisation plane. The squeezed pulses 
are then combined after a tt/2 relative phase shift to form 
a circularly polarised pulse, whose fluctuations in the 
dark Stokes S2-S1 plane are then measured by use of a po- 
larisation rotator and a pair of balanced photo detectors. 
Here we define the Stokes operators as So = N xx + N yy , 
Si = N xx - N yy , S 2 = N xy + N yx , and S3 = iN yx - iN xy 
where Nij = J dz$\(t, z)ifj(t, z) and ^J,^; are the 
quantum field operators for horizontal and vertical po- 
larisations, i = x,y respectively. A half- wave plate is 
used to access different polarisation measurements by ro- 
tating the Stokes operators as Sq = cos (0)Si + sin (0) §2 
in the dark plane [9, 12]. 

For polarisation measurements, the Heisenberg in- 
equality depends on the mean polarisation: A 2 S^A 2 Sj > 

|(5fc)| 2 , where the variance is A 2 S^ = (S^) — (Si) 2 • 
Thus squeezing in the dark plane occurs when A 2 S# < 
|(53)| [9, 12], resulting in the polarisation squeezing met- 
ric M = A 2 Se/\(Ss)\. After linear loss e, the squeezing 
metric becomes M = (1/(S ))[(1 - e)(Si) + e(S )]. 

The simulations are performed by numerically solv- 
ing a stochastic nonlinear Schrodinger equation derived 
in the truncated Wigner representation [13], which has 
proven successful in polarisation squeezing simulations 
for standard fibre [9, 12, 14]. Scaling time and distance 
by the pulse width to and dispersion length zq = t^/fa, 
respectively, and using a propagating frame of reference 
at the group velocity v, we can write the full evolution 
equation as: 
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where we have included third-order dispersion with rel- 
ative strength, B3 = ^3 / 1 ^2 1 ^ an d self-steepening with 
s = The nonlinear response h(r) includes the 

Raman contribution, which has an associated noise term 
r. (See [15] for details.) GAWBS is modelled as Gaus- 
sian fluctuations 77 in the refractive index, with a mag- 
nitude G = 3.2 x 1(T 4 that was determined by a single- 
parameter fit to the low-pulse-energy squeezing experi- 
mental results [11]. 

The field is rescaled as <\> = (vto/n)^ where n = 
I |/(^7^o) is half the average photon number for the 
fundamental soliton at frequency cjo in a medium with 
nonlinear parameter 7. We base the fibre parameters 
on a NKT Photonics NL-PM-750 nonlinear PCF [16] 
at a wavelength of 810 nm. The effective modal area 
of the PCF, 2.0 x 10 -12 m 2 , is much smaller than the 
usual optical fibre, leading to a larger effective nonlinear- 
ity 7 = 91.4 (kmW) -1 . With group velocity dispersion 
/?2 = 12.2ps 2 / km and pulse width to = 0.068 ps, the 
dispersion length is zo = 0.38 m and characteristic fun- 
damental photon number ft = 8.0 x 10 6 . The energy of a 
fundamental soliton is then E s = 2nhuj = 3.92 pj. 

A comparison of the experimental and numerical 
squeezing results, over the range of pulse energies used 
in the experiment, is given in Figure la, lb. Even with- 
out taking into account any linear losses or GAWBS, the 
theoretical squeezing results (crosses) are consistent with 
the experimental data in the overall trends and in the en- 
ergy at which best squeezing occurs (14.6 pJ). In general, 
a larger pulse energy leads to a stronger Kerr effect, un- 
til Raman effects combined with pulse reshaping cause 
the squeezing to deteriorate. Quantitative agreement is 
obtained once we take into account linear losses (dia- 
monds), the finite overlap [11] of the pulses (squares), 
and GAWBS phase noise (circles). The best squeezing 
found numerically was -5.1 dB without phase noise, and 
-3.8 dB with phase noise, consistent with the measured 
-3.9±0.3 dB. 

Significant differences between the simulations and ex- 
periment start to emerge once the pulse energy exceeds 
about twice that needed to generate solitons. In this 
regime, where soliton fission and pulse reshaping lead to 
complex pulse evolution, we can expect that small dis- 
crepancies between simulation and experiment can lead 
to large differences in the squeezing. Moreover, higher- 
order effects may become important, due to pulse narrow- 
ing. Here we considered two such effects: self steepening 
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FIG. 1: (a) Squeezing and (b) anti-squeezing as a function of 
the total pulse energy. Pluses show the 'pure' squeezing, with 
no higher-order dispersive effects (B3 =0, s = 0), GAWBS 
noise (G = 0) or loss. The remaining curves successively add 
in the effects of the 13% linear loss (diamonds), imperfect 
spectral overlap (squares) and GAWBS noise (circles). The 
shading on the latter indicates the effect of a ±0.03 uncer- 
tainty in the experimentally determined spectral overlap [11]. 
Experimentally measured squeezing are given by the trian- 
gles. Vertical dashed line indicates the energy of two soli- 
ton inputs, (c) Squeezing and (d) anti-squeezing for various 
amounts of third-order dispersion B 3 : (pluses), 0.3 (circles), 
0.5 (squares), 1.0 (diamonds), 3.0 (crosses) and 6.0 (stars). 
Other parameters are as in Fig. la-b, with experimental re- 
sults again plotted as triangles. 10 4 stochastic paths were 
used in the simulation to give a maximum sampling error of 
±6.3 x 10~ 2 dB. 



and third-order dispersion. We found the results with 
self-steepening are almost indistinguishable from those 
without. The maximum difference of 0.08 dB is of the 
same order of magnitude as the maximum sampling error 
of ±6.3 x 10 -2 dB in the simulation. 

Since the size of the third-order dispersion is unknown, 
we plot in Fig. lc and Id results for a range of possible 
values, Bs = [0,6]. As expected, it has a stronger ef- 
fect at higher pulse energies, and a value of Bs — 0.3 
may explain some of the faster deterioration of squeezing 
observed experimentally. 

So far, we have only considered the squeezing of a fi- 
bre with a fixed length of L = 1.0 m. Figure 2 shows the 
squeezing optimized over fibre length (up to 30 m). Al- 
though the dependence of fibre length is relatively weak, 
as indicated by the breadth of the shaded area, signifi- 
cantly better squeezing can be obtained by use of longer 
fibres. Overall the best squeezing was found to be -12.6 
dB (without any linear losses or GAWBS noise) using 
soliton inputs and a fibre length of around 5m, compared 
with -9 dB achieved with lm of fibre and higher input 
energies. 

We conclude by noting that, from amongst the vari- 
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FIG. 2: Fibre length (shaded green line, left axis) for optimal 
squeezing with GAWBS noise (shaded blue line, right axis, 
with maximum sampling error of ±0.2 dB) as a function of 
total pulse energy, with (solid line) and without (dashed-dot 
line) GAWBS noise. The green shading shows the range of 
fibre lengths that give squeezing within 5% of the optimum. 
For some of the small pulse energies (between 3.54 pJ and 
4.57 pJ) the trend indicates marginally better squeezing for 
fibre lengths larger than 30 m. Linear loss and finite spectral 
overlap effects have not been included. 10 3 stochastic paths 
were used in the simulation. 

oils physical effects that we have considered, there are 
two that are particularly pertinent for the experimen- 
tal setup considered here: GAWBS noise and the imper- 
fect spectral overlap of the counter-propagating pulses. 
Here we see GAWBS reducing the optimal squeezing 
from -12.6 dB to -8.3 dB (without linear loss), wheres 
for the previous single-pass scheme using conventional 
fibre [9], GAWBS noise had a negligible effect on the 
optimal squeezing. The imperfect spectral overlap intro- 
duces an extra effective linear loss that was not present 
in the previous work. However, if these effects can be 
eliminated in future experiments, we can expect to see a 
gain in the degree of squeezing achieved by use of PCFs, 
over and above what has been observed in standard fi- 
bre. One potential benefit of using PCFs that is yet to 
be explored is the scope for reducing detrimental Raman 
effects, through use of pulses with larger time width. 
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